Stochastic processes with finite correlation time: modeling and 
application to the generalized Langevin equation 

T. Srokowski 

Institute of Nuclear Physics, PL - 31-342 Krakow, Poland 
(February 1, 2008) 



Abstract 



The kangaroo process (KP) is characterized by various forms of the covariance and 
can serve as a useful model of random noises. We discuss properties of that process 
for the exponential, stretched exponential and algebraic (power-law) covariances. 
Then we apply the KP as a model of noise in the generalized Langevin equation 
and simulate solutions by a Monte Carlo method. Some results appear to be incom- 
patible with requirements of the fluctuation-dissipation theorem because probability 
distributions change when the process is inserted into the equation. We demonstrate 
how one can construct a model of noise free of that difficulty This form of the KP 
is especially suitable for physical applications. 

PACS numbers: 02.50.Ey, 05.40.Ca, 02.70.Lq 
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I. INTRODUCTION 



It is frequently assumed that a noise in stochastic equations is not correlated and the 
underlying stochastic process can be regarded as Markovian. This assumption is justified 
only if the time scales involved are large compared to the noise correlation time. On the 
other hand, the noise itself can result from a procedure of fast modes removal. It is well 
known |l]J^] that in such case noise correlations, both in space and time, must arise. For 
some stochastic processes the noise covariance decays fast with time and it can be put in an 
exponential form (the coloured noise ||). Frequently even longer tails of the noise covariance 
are observed. Algebraic covariances appear in the fluid dynamics and linearized hydro- 
dynamics J7j; they are responsible for such phenomena as noise- induced Stark broadening 
H and anomalous nuclear scattering in the framework of the molecular dynamics 0. A 
direct consequence of the algebraic form of the velocity autocorrelation function, falling not 
faster then 1/t, is the infinite value of the diffusion coefficient [Kj. An anomalous diffusion 
process is frequently observed in disordered media where a trapping mechanism leads to 
algebraic distributions of waiting time. For example, some amorphous insulating materials 
(e.g. As2Ses) exhibit a highly dispersive transient photocurrent [|ll|], just due to charge 
hopping between spatially disordered sites. 

Stochastic dynamics driven by a noise different from the white noise obeys the generalized 
Langevin equation (GLE) fl2"yi3f : 



m Mt)_ = _ m f* R{t _ T \ v r T \ dT + F u\ n) 
at Jo 

where F(t) is a stochastic force and m denotes the mass of the Brownian particle. Due to 
the second fluctuation-dissipation theorem (FDT) [|14|,|]1|, the kernel K(t) can be expressed 
in terms of the noise covariance Cp = (F(0)F(t)): K(t) = Cf/ttiT, with the temperature T. 
The Eq.(|l|) can be handled as a usual Volterra equation ||16|| . Assuming the initial condition 



v (0) = 0, the general solution can be expressed in the form of a stochastic integral [17 



v(t) = m- 1 f R(t - r) F(t) dr, (2) 
Jo 
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where the Laplace transform of the resolvent R(t) is given by the equation 

R(s) = l/[s + K(s)]. (3) 

From the Eq.(Q), expressions for some average quantities follow. For the velocity variance, 
we have 

( v 2 ) s (t) = m' 2 /* /* R(t - r) R(t - r') C F (\r -r'\) drdr', (4) 
J o J o 

where the average ( )g is taken over an equilibrium ensemble with some stationary proba- 
bility distribution. The FDT ensures that asymptotically, for large time, the system reaches 
the equilibrium value (v 2 )s = T/m (the equipartition energy rule). For the velocity auto- 
correlation function, in turn, we have simply: 

C v (t) = (v(0)v(t)) s =T/mR(t). (5) 

The assumption about a form of noise covariance is sufficient to calculate some average 
quantities. In order simulate stochastic trajectories from the Eq.(0), one needs a concrete 
physical process which could serve as a model of the noise. For example, for the exponential 
covariance it could be the well-known Ornstein-Uhlenbeck process, which, due to the Doob 



theorem [18], is very important if amplitude distributions are Gaussian. A broad class of 
stochastic processes known as "kangaroo processes" (KP) fl]| is especially interesting. One 
can construct the KP for an arbitrary, given form of covariance. The KP is particularly well 
suited for problems involving algebraic, scale-invariant dependences. Long tails of constant 
value of the process in the step-wise structure of the KP make possible to preserve the 
memory about this value for a sufficiently long time to produce such slowly decaying form 
of the covariance. Due to that structure, the KP resembles stochastic, dispersive transport 
processes in disordered media, e.g. the hopping time distribution [f20|-|22f. In the framework 
of random walk processes, a pattern of long straight-line segments is typical for Levy flights 

EMI- 

This paper deals with random noises possessing various covariances: exponential, 
stretched exponential and algebraic, and expresses them in terms of the KP. The most 



important properties of KP are summarized in Sec. II; we also derive there formulas refer- 
ring to those forms of the covariance. In Sec. Ill we consider the application of the KP 
as a model of random force in the GLE. The most important results are summarized and 
discussed in Sec. IV. 



II. THE KANGAROO PROCESS 



The kangaroo process |T9| is a step-wise, discontinues random function. The value of the 
process, m(t), is determined at subsequent random jumping times ti, t 2 , ■ ■ ■■ The jumping 
frequency v(m) depends on the value of the process itself and m remains constant between 
jumps. We introduce also the interval length as a reciprocal of the frequency: s = \ jv. Due 
to some physical applications, this quantity can also be called "a free path". The KP is a 
stationary Markov process and can be defined by the probability density p(m, t) satisfying 
the following Fokker-Planck equation 

^-p(m,t) = v{m) \-p(m,t) + - jrjf-^r n , , / v(m')p(m' , t)dw! J (6) 
at V J v{m')PKp{jn )dm' J J 

where Pjcp(m) denotes a stationary probability distribution of m(t). The interval length 
s is also a stochastic quantity. Its probability distribution P(s) is connected with P KP {m) 
by the relation P(s)ds = 2PKp(\m\)d\m\. We assume that Pxp{m) and u(m) are even 
functions of m. This assumption allows us to get a simple expression for the covariance 
of KP C(t) = (m(t)m(0))s where the average is taken over the stationary probability 
distribution P KP {m): 

/+oo 
m 2 P KP (m)exp(-u(m)\t\)dm. (7) 
-oo 

We want to derive an expression for u(m) for a given covariance C(t) and an amplitude 
distribution P^-p(m). Let us assume that vim) is a monotonic increasing function of |m| 
and v[oo) = oo. Then we can change the integration variable in Eq.(^) and obtain the 
Laplace integral 
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C(t) = 2 m 2 P KP {m) — exp{-v\t\)dv. (8) 
Jj/(0) az^ 

Therefore v(m) is a solution of the following differential equation 

^- = 2m*P KP {m)/C{v) (9) 

where C(V) denotes the inverted Laplace transform of C(t). 

Solution of the Eq.fl^) allows us, in principle, to generate a stochastic time series of the 
process with a given, quite arbitrary covariance and with an arbitrary distribution Pxpim). 



A. Exponential covariance 

The KP for the exponential covariance 

C{t) =z/ exp(-z/ot), (10) 

where vq =const is a reciprocal of the correlation time, is called the Kubo- Anderson process 
|26f . The jumping times are uniformly distributed in the interval (—00,00) with a m- 



independent density Vq, according to the Poissonian distribution. Therefore the intervals of 
constant m are distributed exponentially: P(s) = u exp(— u s). The value of the process 
m(t) may be chosen according to an arbitrary distribution Pxpim). In that sense, the 
distributions of m and s are independent of each other. That property holds only for the 
Kubo- Anderson process; for a general KP m and s are interdependent. The Fokker-Planck 
equation @ takes a simpler form for the Kubo- Anderson process: 

J^o(m, t) = u (-p(m, t) + P K p(m) J p(m', t)drrA . (11) 

Some forms of P^p(m) are distinguished. The simplest choice is P^pim) ~ 5(m — u ) + 
S(m + vq) and corresponds to the dichotomous noise (the random telegraph process) |27| . 
Due to the central limit theorem, the Gaussian distribution of Pkp(itl) is of special impor- 
tance. The Kubo- Anderson process with that distribution resembles the Ornstein-Uhlenbeck 
process. However, both processes are not identical [p8 ]; the Fokker-Planck equation for 



the Ornstein-Uhlenbeck process, j^p(m } t) = (m + Dv -£^ p(m,t), differs from the 

Eq.©. 

B. Stretched exponential covariance 

The exponential distribution of interval lengths P(s), a distinctive feature of the Kubo- 
Anderson process, can also characterize stochastic processes with the covariance form other 
than exponential. Let us consider the covariance possessing the shape known as the 
"stretched exponential". This function is the following: 

C(t) = exp(-c^) (12) 

where a =const and < 7 < 1. In the present study we assume 7 = 1/2. The stretched 



exponential function describes relaxation phenomena in random systems |29,30|| and can 



be attributed to a dispersive transport of mobile defects in the glass [[JTJ. The dispersive 
transport is characterized by the infinite average time between subsequent hops. Moreover, 
the velocity autocorrelation function of a particle inside the Sinai billiard with the finite 
horizon is also given by the Eq.(O) |R3|. 

We want to find expressions, useful in practical applications, allowing us to generate time 
series possessing the required covariance. The inverted Laplace transform of fll2D reads 



^) = ^^ 3/2 exp(-a 2 /(4^)). (13) 

2 \J IT 

In order to simplify the differential equation @, we take the amplitude distribution in the 
form 

2 

Pft-p(m) = — |m|~ 5 exp(— TT^m" 4 ), (14) 

7T 

where m E (—00,00). The distribution possesses maxima at m — m max = ±[4/(57r)] 1 / 4 
and it is very small near m = 0. Inserting Eq.fllBD and Eq.flTj]) into the Eq.(^), we get the 
differential equation for vim) in the form 
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= 8 /( v ^a)m- 3 z/ 3 / 2 exp[-^(z/- 1 - 4/(a 2 vrm 4 ))]. (15) 

One can easily check that Eg. (|T5|) is satisfied by the function 

u(m) = l/s = (7ra 2 /4) m 4 (16) 

and the interval lengths distribution is indeed exponential: 

P( S ) = ^exp(-^ S ). (17) 

The direct relation between s and m follows from the Eq . (|i"6|) . From that equation we 
conclude that s G (0, oo) and long intervals correspond to the values of m close to zero; such 
events are extremely rare. Technically, a time series m(t) can be constructed by sampling 
subsequent intervals s from the distribution P(s). Then corresponding process values are 
evaluated by means of the Eq. (|T6|) , taking into account, in addition, that both signs of m 
are equally probable. 



C. Algebraic covariance 

Let us now consider the KP possessing the power-law covariance that we express in the 
following form 

<7(t) = r( 7 ) *-> ( 7 >0). (18) 

The jumping frequency vim) one can derive from Eq.(||), similarly as for the stretched 
exponential case: 

l/ 7 

(19) 

We assume the amplitude distribution in the algebraic form: Pxpivn) ~ vrC a (a = const), 
modified in order to satisfy the normalization condition 2 J °° PKp(m)dm = 1. Generally, 
two different forms of Pxp(^) are possible. 

In the first case we cut off the large values of |m|: 

7 
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\m\ 



m 



l P KP (m')dm' 



hr^(\m\/a) a \m\ < a 
P KP {m) = { 2a ~ (20) 

|m| > a 

where a =const is an additional parameter. Due to the condition u(oo) = oo, a must be a 

large number; finiteness of a results in a deviation of the covariance from the assumed form 

(|18|) near t = and in removal of the singularity. Moreover, the normalization condition 

implies a < 1. Inserting P^-p(m) to the Eq. (|l9"D gives us the expression for the relation 

between m and s: 



s = 

V 



- = ( ia ')- lh a {1 - a) ^\m\-^- a ^9(s - e) (21) 



where we have introduced a constant a' = (1 — a)/ (3 — a). Finiteness of the parameter a 
imposes a restriction on the lower bound of the interval length e: s G (e, oo). The smallest 
interval length e is related to that parameter by: e = (0/7) _1 / 7 a~ 2 ^ . The probability 
distribution of interval lengths takes the algebraic form 

P( s ) = {ja'f-^a-^'s-^ 1 - 1 . (22) 

From the Eq. (p2|) some restrictions on possible asymptotic behaviour of P{s) follow. The 
slowest decay rate for large s occurs for a close to 1: P{s) ~ s . On the other hand, the 
distribution P(s) falls the most rapidly, as ~ s~ 7_1 , for a — > —00. 

The other possibility is to cut off the distribution Pxpim) near m = 0: 

!0 \m\ < a 

(23) 
^-(\m\/a) a \m\ > a 

where a > 1. One can show that the minimal interval length is finite (non-zero) if a > 3. 
The interval lengths distribution can be obtained similarly as for the case (p0|). The result 
is slightly more complicated: 

-2/(3-a) 

Asymptotically the distribution approaches s~ 7_1 , independently of a 



P( s ) = [ — * a i-« s -i + a 3 ~ a ) s-^-\ (24) 
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III. APPLICATION TO THE GLE 



We will now consider the Eq.(|I|) for which the random force F(t) is modeled by means of 
a concrete stochastic process, possessing a given covariance, and simulate stochastic trajec- 
tories of the Brownian particle by means of a Monte Carlo method. Therefore, two equations 
are to be solved simultaneously: the original GLE and the second one, describing an ad- 
joined process, in the form of some KP. Accordingly, in the following we assume F(t) = m(t). 



A similar approach used to be applied |x|o3J] to deal with stochastic equations driven by 
the coloured noise; the adjoined process constitutes in that case the Ornstein-Uhlenbeck 
process. Independently of those methods, which model the noise as a stochastic, random 
process, one can introduce some deterministic system possessing a required autocorrelation 
function. Value of the noise at a given time is then determined by the evolution of dynamical 
equations of motion. In this way Shimizu |35| solved both ordinary Langevin equation and 
GLE, representing the noise by a chaotic map. Similarly, the Sinai billiard has been used to 
model a noise with long-time covariance in the Langevin equation |36| . Certainly, the first 
two moments do not determine the noise uniquely and the choice of a model involves addi- 
tional assumptions, especially for non-Gaussian and not exponentially correlated processes. 
The form of the noise must be decided according to physical requirements of a concrete 
application. Nevertheless, some important quantities, like the autocorrelation function of 
the Brownian particle velocity (||]), apparently depend only on the noise covariance. 

However, for some stochastic processes description of those quantities provided by the 
equations collected in the Introduction does not agree with the simulation results. Solving 
the GLE to obtain the Brownian particle velocity requires the value of the noise F(t), 
determined by the adjoined process, at a given time. It has been demonstrated recently 



|37| that this requirement modifies probability distributions - the adjoined process looks 



differently when inserted into the equation. Consequences of that change for the Monte 
Carlo simulation results may be important. In particular, one can expect a violation of the 
FDT which manifests itself in an improper asymptotic behaviour of the velocity variance 
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{v 2 (t))s ■ Determining of that quantity can indicate whether the expected equilibrium state 
is reached and the equipartition energy rule holds. 

The Brownian particle velocity is to be determined from the stochastic integral (Q) in 
which step-wise, constant values of F(t), Ff~, follow from the length of subsequent time 
intervals s^: 

v(t) = m~ 

where tk — tk-i = Sfc and to = 0. Sampling of intervals s& continues as long as the time t 
is reached: J2i Sk < t and Sk > t, where n is an arbitrary integer. The interval length 
distribution P(s) is a natural quantity of interest. The last, n+1, interval corresponds to the 
time t [p8[| . It is clear that the distribution of lengths of that interval, the "effective" interval 
distribution P(s,t), cannot be identical with P(s); a simple consideration reveals e.g. the 
enhanced probability of choosing longer intervals. Generally, that modified distribution can 
depend on t. One can express it in terms of the cumulative distribution function $(s, t) p7| : 

P(s,t) = ^(s,t), (26) 

where 

(J*_.S(x)dxfl_P(£)dt for 0<s<t 

$( S ,t)= (27) 
[SiS{x)dx!^P(Z)dt + !;P{Z)dZ for 8 >t 

The form of auxiliary function S(x) follows from the normalization condition 

rt roo roo 

S(x)dx P(m+ P{Z)d£=l. (28) 

JO Jt-x Jt 

For the Kubo- Anderson process, the modification of the interval distribution P(s), ex- 
ponential in that case, is of minor importance because Pkp{F) is independent of P(s) and 
the interval lengths do not influence process values. Consequently, results of the simula- 
tions agree with general predictions implied by the FDT. Other forms of noise covariance 
require taking into account the modified distribution P(s,t). The case of stretched expo- 
nential covariance provides a simple but non-trivial example. It involves also the Poissonian 
distribution of jumping times. 

10 



n 

E 



tk- 



R{t -t)cLt + F n+1 / R(t - t) dr 



(25) 



A. Stretched exponential noise covariance 



We assume the noise covariance in the following form: 

C F (t) = 2mT/a 2 exp(-aVi). (29) 

The Eqs.(Q) and fl25|) express the solution of GLE in terms of the resolvent R(t). The Laplace 
transform of that function, R(s), is given by the Eq.([3D and kernel has the form: 

where erf(x) denotes the error function. To obtain the resolvent R(t) we need to evaluate 

the inverse Laplace transform from R(s): R(t) = \ f^^t* R(z) e tz dz. The integrand 

possesses two conjugate simple poles and a cut along the negative real axis. Evaluation of 

the contour integral produces the following result 

_ at 4 r°° x 2 exp[x 2 - a 2 t/(Ax 2 )] dx 

R{t)=e [c lS mbt + c 2 cosbt)--j=J o [(2x 2 + a 2/( 4x 2)) exp(x 2) _ 2 0Fx 3 erfi(a;)] 2 + 4ttx 6 ; 

(31) 

the imaginary error function erfi (x) = — % erf (ix) can easily be calculated by the following 
expansion: 

2 oo x 2n+1 

The constants a and b denote the real and imaginary parts of the pole of R(z), respectively: 
z = — a — \b\i; they have to be evaluated numerically. c\ and c 2 can be found by the standard 
residues analysis. For a = 1 the constants are the following: a = 0.207094, b = 0.440963, 
C\ = —0.127752, and c 2 = 0.593952 . Fig.l presents the function R(t) for a = 1 and a = 2. 

We wish to perform the Monte Carlo simulation using the noise defined by the process 
(fnj). According to (|T7j), the interval distribution is exponential: P(s) = /3exp(— (3s), where 
(3 = a 2 /4. We expect that taking into account of the modified form of the distribution may 
be important for the simulation results because m and s are connected. That distribution 
can easily be found in this case. From the Eq . (|28D we obtain S(x) = /3; finally we get: 
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f /? 2 sexp(-/3s) for 0<s<t 

P(s,t) = \ (32) 
{(3(1 + f3t)exp(-ps) for s>t. 

The function (^) is presented in Fig. 2. The distribution is discontinues. It depends on 

t but this dependence dwindles exponentially with time; the left-hand branch (s < t) is 

time-independent. Therefore, asymptotically the process becomes stationary. Nevertheless, 

P(s, t) possesses the mean value twice of that for P(s) [^] and simulation results must reflect 

that. Indeed, the velocity variance obtained from the Eq. (|25|) , shown in Fig.3, stabilizes at 

a lower value then that predicted by the equipartition energy rule (T/m). It is so because 

long intervals correspond to small values of the noise amplitude, according to the Eq. (|T6|) . 

For short times, in turn, the result of the simulation agrees with the general prediction (£|), 

also shown in the Fig.3, because then the branch s > t dominates the distribution (|32l) and 

the dependence on t is weak. 

It is possible to construct some KP which does not change when inserted into the GLE 

and which produces the equilibrium state in agreement with the FDT. For that purpose we 

single out some subset of kangaroo processes, a restricted KP (RKP) , defined in the following 

way. We choose subsequent intervals Sk in the step-wise evolution of KP, according to the 

distribution P(s), and assume that t corresponds to n+1 interval, i.e. S n = S\+S2 + - ■ -+s n < 

t and S n+ i > t. Let S t = S n+ \ — t > 0. We call some KP "restricted" if St < d for a given 

d; n is an arbitrary integer: n G (0, oo). The probability distribution of the last interval 

lengths we denote by P*(s,t;d). Obviously, P*(s, t;oo) = P(s,t). On the other hand, in 

the limit d — >• the last interval lengths obey the original, time-independent distribution 

P{s) but only for s < t because longer intervals are excluded by construction. Therefore 

the distribution P(s) can be recovered at large times if those intervals are negligible. We 

get then the following theorem: If the probability that intervals in a sequence Sk are larger 

than t asymptotically vanishes, then 

lim limP*(s,t;tZ) = P(s). (33) 

t^co d— >0 

For the distribution (|32|) the required probability vanishes with time and the theorem can 
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be applied. In practice it can be done easily by choosing some small d and sampling intervals 
from the distribution P(s), in the same way as before, as long as the time t is reached. 
Then all sequences of intervals for which St > d are rejected. Fig.3 presents the result of 
such calculation for d = 0.01. At short times, the velocity variance differs substantially 
from the other results shown in the Figure because the branch s > t is then essential but 
asymptotically it approaches the value T/m, in accordance with the equipartition energy 
rule. 



B. Algebraic noise covariance 

Finally, we consider a power-law form of the covariance: 

C F (t) ~ r 7 (0 < 7 < 3). (34) 

The KP we apply to model the noise is defined by Eq. (p0|) ; we implement the simplest case 
a = 0. Then the amplitude distribution is a constant (except very large \F\): Pkp(F) = 



7e 7 /3/2 where e is the smallest interval length. The interval distribution follows from 
Eq.©: 

P( s ) = 2 e T/3 s -i-7/3 0( a _ e ). (35) 
3 



The resolvent R(t) for the noise with covariance fl34j) can be evaluated by means of 
similar methods as for the stretched exponential covariance. Results can be found in Ref. 
40fl for 7 = 3/2 and in Refs. [37 ,fflH for 7 = 1. Now we want to calculate the modified 



distribution P(s,t). First, we have to solve the Eq . fl28|) which assumes the form of Abel's 
integral equation 

* S(x)(t - x)-~ <ri dx + r 7/3 = e~ 7/3 . (36) 



The solution reads 



f -7/3 

g W = r(i- 7 /3)r (7 /3) * 7/ ~ S{X) - (37) 
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After evaluation integrals, we obtain from Eqs. ( fffi) and (p7|) the expression for the required 
distribution: 

S-7/3-1 



P(s,t) = { 



r(i - 7 /3)r( 7 /3) 

t 7/3 



£ 7/3 _ ( t _ s )7/3 
7/3 



for e < s < t 



(38) 



+ 



7e' 



-7/3-1 for 



S > t. 



r(i- 7 /3)r( 7 /3) ' 3 

The distribution P(s, t) for 7 = 1 is presented in Fig.4. The picture is markedly different 
from that obtained for the exponential case (Fig.2); the right branch, corresponding to the 
intervals s > t, does not vanish with time but gets larger comparing to the left branch. 
The entire distribution shifts with time towards large intervals because the average interval 
length is infinite. 

Therefore the effective distribution P(s,t) differs substantially from the P{s), and the 
stochastic process generated by it must possess different properties. First, let us recalculate 
the covariance that, in general, can depend on an initial time t : C(t, t ) = (F(t )F(t +t)) 5 , 
where the process F(t) is to be determined by the simulation. Technically that means that 
for a given to, one produces a sequence of intervals to reach the time to+t. Then one evaluates 
corresponding values of the process F(t) using Eq.(^). Then the process is governed by the 
distribution P(s,t). The expression for the covariance follows from the Eq.(^): 



C(t,t ) 



-27/3 



exp(— tj s)P(s, t )ds. 



(39) 



Evaluation of the integral gives the following result 
C(t,t 



3^7/3^-1/3 

[f /3 t-^h,t/e) 



r(i-7/3)r(7/3)' 



t l/2- 7 /6 t( _ 1 _ 7)/2exp( _ t/2to)r(1 + 7/3) ^ 7/6 _ 1/W2 ( ^ o)]) 



(40) 



where 7(0, x) denotes the incomplete gamma function and W a $ (x) stands for the Whittaker 
function |4"2| . Certainly, the above result is different from our starting covariance (0); the 
most striking feature of the function C(t, to) is the dependence on to which does not diminish 
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with t. The variance of the process, a 2 (to), can be found by inserting t = into the Eq. (fl0|) . 
Let us consider two examples. The case 7 = 1 has been discussed in Ref. |57H ; the final 
expression for the variance is the following: 

° 2 (t ) = r(1/ e 3) ^ 2/3) [3 In 3/2 + vrv/3/6 + ln(t /e)] t^ 3 (t » e). (41) 

Algebraically correlated stochastic processes for 7 = 3/2 are especially important. They 
have been extensively studied in connection with the Brownian motion in a viscous fluid 
4l||4]-|6|,[40| . In this case the variance of the process reads 



- 2 (to) = (^) 1/3 (^to 1/2 + t 3/2 ). (42) 

Therefore, in contrast to the original variance (calculated with the distribution P{s)) a 2 = 
7~ 1 e~ 7 , the effective variance a 2 is time-dependent and tends to zero. That behaviour is a 
direct consequence of nonstationarity, i.e. of the time-dependence of the distribution P(s, t). 
The decline of a 2 (to) means that the effective temperature of the system drops to zero when 
we insert the process into the stochastic integral (EH) as a model of the noise. Consequently, 
the Brownian particle velocity variance also must dwindle with time. Indeed, a direct Monte 
Carlo simulation confirms this conclusion, as it has been demonstrated in Ref. [37J for 7=1. 

The RKP can be constructed also for algebraic correlations. However, the theorem (|33"D 
cannot be applied because of the strong time-dependence of the distribution P(s,t). The 
probability that the interval length is larger then t does not decline with time: P(s > t,t) = 
J t °° P(s, t)ds = 3/[ 7 r(l - 7 /3)r( 7 /3)] =const. 



IV. SUMMARY AND DISCUSSION 

The kangaroo processes represent a broad class of random functions characterized by 
various forms of the covariance. Therefore, they provide an opportunity to model physical 
stochastic processes possessing an arbitrary covariance and, in addition, a quite general 
amplitude distribution. The KP is step-wise - the value of the process changes according to 
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some jumping frequency which, in turn, depends on that value. Some physical phenomena 
exhibit a similar, step-wise behaviour and the KP is a natural process to model them. An 
important quantity is the distribution of intervals of constant process value P(s), uniquely 
connected to the amplitude distribution. In this paper, we have discussed two forms of 
that distribution: the exponential and the algebraic ones. We have demonstrated how one 
can generate algebraically correlated processes using the KP with some algebraic P(s). The 
exponential form, in turn, is suitable to represent both exponential and stretched exponential 
correlations. 

The GLE has been solved using KP as a model of the noise. The problem has been 
considered as a juxtaposition of two random processes: of velocity of the Brownian particle, 
described by the GLE, and of the adjoined KP. In the framework of that approach, the 
force F(t) in the stochastic integral (0) is determined by looking for a value of the KP, 
independently evolved, at a given time. Such procedure changes probability distributions of 
the KP and Monte Carlo simulated solutions of the GLE are not in agreement with results 
predicted by general analysis, founded on the FDT. Results obtained from effective prob- 
ability distributions, among which the interval distribution P(s,t) is the most important, 
does not correspond to the equilibrium state consistent with the equipartition energy rule. 
For algebraically correlated processes, even nonstationarity effects emerge. Are those results 
involving effective distributions a necessary consequence of modeling of the noise by means 
of adjoined random process in the form of KP? If the stochastic force we want to insert 
into the stochastic integral (|2|) represents some independent physical process, the modified 
probability distribution (|26|) has to be taken into account. In such case we must expect 
an apparent violation of the FDT, despite proper definition of the kernel. Reversely, any 
information about properties of the noise, extracted from GLE solutions, always refers to 
the effective appearance of that process in the GLE. 

On the other hand, for some cases there is a possibility to construct the model in such 
a way to avoid any modification of the distributions and to preserve consistency with the 
FDT. For that purpose, one can use some specific version of the ordinary KP - the RKP. 
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The idea is very simple: if some interval in the step-wise evolution of the KP ends exactly at 
a point corresponding to the required time, nothing has to be modified. We can imagine the 
KP as a "clock" with a variable frequency given by the distribution P(s). In that picture, 
applying of the RKP with d — means a synchronization of that clock in respect to the 
physical time in the GLE. An important limitation of the synchronization procedure consists 
in the fact that interval lengths of the RKP are always finite, not larger than the time t, 
and P(s) usually possesses infinite tails. If, however, those tails decline sufficiently fast, 
intervals larger then t become negligible. It has been demonstrated that the exponential 
interval distribution possesses this property and a Monte Carlo algorithm utilizing the RKP 
can easily be constructed - simulation results indeed correspond to the equilibrium state 
predicted by the FDT. The KP in this form is especially suitable for physical applications. 
We note, however, that the RKP with d — is not completely independent of the GLE: the 
synchronization introduces a coupling. The case of stretched exponential covariance is only 
the simplest non-trivial example of application of KP with the Poissonian interval lengths 
distribution. A straightforward generalization, allowing for the other power-law dependences 
m(s), produces KP's with covariances given by the Bessel functions K v (y/i). 

One cannot expect that every shape of the covariance function may be modeled by some 
sufficiently steep form of P(s) and, therefore, the RKP is always a proper tool. Processes 
possessing covariances with long tails are characterized by long intervals (long free paths). 
Indeed, in Sec. II we have indicated a strong limitation of admissible shapes of P(s). The 
algebraic tails of P(s) are essential and very long intervals given by the effective distribution 
P(s,t) are not negligible for large times; the distribution itself remains time-dependent. 
For that reason, the RKP which consist in cut off of long intervals, cannot work as a 
model of algebraic covariances. In the other words: due to divergence of moments of both 
distributions, P(s) and P(s,t), the branch corresponding to the intervals s > t is important 
for large t and the synchronization cannot be achieved. 

Necessity of taking into account the modified form of probability distributions is not 
restricted to the GLE but relates to any stochastic equation, e.g. the ordinary Langevin 
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equation, if the random force is modeled by the adjoined KP. Monte Carlo simulations can 
be useful for some generalizations of the diffusion equation, in particular for the Burgers 
equation ||44j| . That nonlinear equation, possessing a broad spectrum of applications in the 
fluid dynamics, is characterized by long range noise correlations both in space and time 
(TJ. Also that noise can possibly be modeled by the KP. 
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Figure captions 

FIG. 1. The resolvent R(t) fl3~l"|) for the stretched exponential shape of noise covariance 
as a function of time. 

FIG. 2. The time evolution of the effective interval distribution P(s, t) ([32]) corresponding 
to the exponential form of the original distribution P(s) with (3 = 1. 

FIG. 3. The velocity variance obtained from the Monte Carlo simulation of GLE solutions 
(p5|) for stretched exponential form of noise covariance. The noise has been modeled by 
the KP (solid line) and the RKP with d = 0.01, using theorem (|33|) (dashed line). 1000 
trajectories have been calculated for either case. The variance calculated from Eq.(|4]) is also 
presented (dots). The temperature T = 1, the particle mass m = 1, and the parameter 
a = 1. 

FIG. 4. The time evolution of the effective interval distribution P(s, t) fl3~8|) corresponding 
to the original distribution P(s) in the form (|35|) with 7 = 1. The minimal interval length 
e = 0.01. 
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